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I. INTRODUCTION 

The Bethe-Salpeter (BS) formalism is a venerable method for obtaining information on bound state systems in a 
covariant framework lj. It is, however, technically challenging and most applications to hadronic physics have been 
made in simple models [2j or in two-body systems [3], Only recently has the rather formidable extension to nucleons 
been undertaken 0]. 

Because the Bethe-Salpeter formalism is covariant, it is especially useful for describing dynamical quantities such 
as form factors and distribution functions. It is also a many-body approach and thus can incorporate chiral symmetry 
breaking, which is crucial to obtaining reliable predictions in the light hadron spectrum[3]. 

Unfortunately, these benefits come at a cost. The Bethe-Salpeter formalism hinges on the assumed form of a two- 
body irreducible scattering kernel, and this form must be obtained with some truncation - a truncation that must be 
made with no apparent small parameter in sight. This theoretical issue can become acute because truncated kernels 
often lead to bizarre predictions such as 'ghost states' [2], nonsensical spectra, or unphysical parameter-dependence of 
observables. Even this undesirable situation can be a distant goal. Often one simply postulates a kernel and tests 
the accuracy of its predictions. Of course this removes all putative connections to Quantum Chromodynamics (QCD) 
and the Bethe-Salpeter formalism devolves to a (quite sophisticated) quark model. 

Here we explore the extension of the Bethe-Salpeter approach to the gluonic sector of the strong interactions - 
namely the glueball spectrum. To our knowledge this is the first attempt at computing glueball masses with the 
Bethe-Salpeter formalism. It is also a rare example of a hadronic computation that uses QCD to build the interaction 
kernel [5]. Thus we seek to leave the realm of pure model-building and enter one in which field-theoretic truncations 
become the relevant issue. We will, for example, assume that dressed versions of the lowest order kernel dominates 
the interaction. Thus the ingredients necessary for this investigation are gluon, ghost, and quark propagators and the 
three-gluon, four-gluon, quark-gluon, and ghost-gluon vertices. It is clear that the problem is formidable, even at this 
simple level of truncation. 

Ensuring the fidelity of the truncations employed is not straightforward. We are therefore fortunate that a variety 
of lattice gauge computations of two-point functions have been performed recently. Our strategy will be to model 
vertices, employ these vertices in Schwinger-Dyson equations for two-point functions, verify the solutions against 
lattice computations, and then use the resulting suite of n-point functions in the glueball Bethe-Salpeter equation. 
Unfortunately, this is difficult to implement in practice because the gluon two-point function is notoriously truncation- 
dependent [6l [7] . We therefore choose to use the lattice ghost and gluon propagators directly in the Bethe-Salpeter 
equation. 

We start by defining the various propagators and vertices used in this work. The two-body Bethe-Salpeter equation 
for gluonic states is then developed and applied to even charge conjugation scalar glueballs. Model parameters are 
tuned against lattice gauge theory results and finally, we test the fidelity of the computation by comparing model 
computations of the relevant two-point functions to lattice results. Although a reasonably good spectrum is obtained, 
the comparison to two-point functions indicates that substantial truncation error occurs and further work is required 
to obtain robust QCD-based predictions of hadronic properties in the Bethe-Salpeter formalism. 

II. PROPAGATORS AND VERTICES 

Quark, gluon, and ghost propagators are defined in terms of dressing functions as follows: 
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The scattering kernel that underpins the Bcthc-Salpeter formalism can be formally written in terms of Schwinger- 
Dyson equations. Of course these equations must be truncated at some point. It is traditional to organise the 
equations in terms n-point functions and truncate at the two- or three-point level. We will follow tradition here 
and thus construct the two particle irreducible four-point function via dressed tree-level interactions. Details will be 
provided in the next section - here we focus on the dressed vertices. 



A. Model Vertices 



In the case of abelian gauge theory, the Ward-Takahaski identity can be used to completely constrain the longitudinal 
dependence of the full fermion-gauge boson vertex. The result is given in terms of the fermion dressing functions, A 
and B, and is called the Ball-Chiu vertex[S]. The situation is more complicated in the nonabelian case because the 
ghost dressing function and the quark-ghost scattering kernel, K g appear in the Slavnov- Taylor identity: 

p£r£(pi,Pa,P3) = ig T a h(p 3 )[S- 1 (-pi)K g (p 1 ,p 2 ,p 3 ) - K g (p 2 ,p 1 ,p 3 )S- 1 (p 2 )]. (4) 

The resulting vertex has a rich structure that was studied by Aguilar and Papavasilliou[9]. They found that an 
adequate approximation to the full longitudinal structure of the vertex can be obtained by neglecting the quark- 
gluon scattering kernel and multiplying the Ball-Chiu vertex by the ghost dressing function. We make the further 
simplification of retaining only the 'central' term (proportion to a gamma matrix) of the Ball-Chiu vertex: 

r>i,p 2 ,P 3 ) = -igoT a 7 M h n (pi - pa) \[A{px) + A(p 2 )}. (5) 

A well-known shortcoming of the Ball-Chiu approach is that the transverse part of the vertex function remains 
unspecified. This has been partly remedied in the Abelian case by requiring multiplicative renormalisability in the 
asymptotic regime [TQ]. No such Ansatz exists for the nonabelian case and we choose to neglect this contribution 
henceforth. However, Aguilar and Papavasilliou note that the correct perturbative asymptotic behaviour for the 
quark propagator is not obtained with this approximation - presumably due to the missing transverse components of 
the vertex. They then observe that supplying an extra power of the ghost propagator in Eq. [5](n = 2) allows recovery 
of the expected asymptotic behaviour. We therefore choose to allow for varying values of n in the following. 

The ghost-gluon vertex has been considered by von Smekal et al.\6\ who neglect the ghost-ghost scattering contri- 
bution to the vertex and retain only the reducible part of the four-ghost scattering matrix element. The result is an 
Ansatz of the form 



K bc (q,p) = g r bc 



h{k) 



U(p) ) 



(6) 



where q is the incoming ghost momentum, p is the outgoing ghost momentum, and k = p — q is the incoming gluon 
momentum. 

The authors of Ref. [5] also examine the three-gluon vertex. They employ the ghost scattering kernel implicit in 
Eq. [6] to obtain an Ansatz as follows: 



r(123) = 5o / 123 
where 
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This form ignores the unconstrained transverse components, but these are known to be suppressed [TT]. Since one 
expects A_ <C A + we simplify the vertex further be neglecting the terms proportional to A- . 
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An alternative model is advocated by Pennington and Wilson [7 who employ an approximate solution to the 
Slavnov- Taylor identity with a bare ghost-gluon scattering kernel: 



(9) 



Vertex models for the four-gluon vertex have been rarely considered [12]; we simply use the perturbative vertex in 
the following. 

B. Gap Equations 

The propagator dressing functions of Eqs. [3] are obtained by solving the two-point Schwinger-Dyson equations. 
These require renormalisation which we implement with Z-factors as defined in the following Lagrangian: 
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FIG. 1: Quark and Ghost Gap Equations. 

The exact quark and ghost gap equations are shown in Fig. [T] As discussed above, these are truncated by 
employing vertex models. Projecting the full quark propagator onto its two components, working in Landau gauge, 
and performing the Wick rotation to Euclidean space yields the coupled gap equations: 



7>, 2 , -/ ., 2 7 2 7 r , I d A q F B(q 2 E )G E (Q E ) 

B(p E ) = ZpnriQ + 3g Z F Z A C F / ttt-tjW - 



A 2 q 2 E + B 2 



(11) 



A(p|) =Z F - g 2 Z 2 Z A C F 



d*q E w A(q E )G E (Q E ) f2q%(l-x 
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A 2 q 2 + B 2 



Pe J 



(12) 



and 



HpD 



= Z c - g 2 Z 2 Z A C A J ^l-x 2 dxVG E {Q E )h{q 2 E )- 



(13) 



Here Q = p — q and the Euclidean nature of these expressions has been emphasised with the subscript, which we 
subsequently suppress. By convention G E (Q E ) = —G(Q Qe), and B(p E ) = B(p — >• p E ). Factors due to vertex 
models are given as 



h(Q) h(Q) 
h(q) + h(p) 



(14) 



and 



W = h n (Q)~(A(p)+A(q)). 



(15) 
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Rainbow-ladder results are recovered when h — > 1 and A — > 1 in these expressions. 

The gluon gap equation, shown in Fig. [2j has a rich structure that makes obtaining definitive numerical solutions 
notoriously difficult [6j [3 [9] . We therefore adopt the simple expedient of employing a fit to the lattice gauge results. 
It is important to note that this is not just a simple matter of finding a parameterization that mimics the data in a 
reasonable way because solving the glueball Bcthe-Salpeter equation requires knowledge of the gluon propagator in 
the Minkowski regime. This will be discussed in greater detail below. 
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FIG. 2: Gluon Propagator Schwinger-Dyson Equation, 
full vertex. 



Solid circles represent full propagators. The open circles represent the 



We chose to renormalise the quark and ghost gap equations by setting 



g Z F Z 



1/2 _ 



= 9- 



However, in Landau gauge 



5o Z C Z 



1/2 



(16) 



(17) 



holds to all orders[13 . Thus gauge covariance implies that Z c (/i, A) = Zp(/i, A), where A is the scale associated with 
the regularisation scheme. Truncation and model error break gauge covariance and thus the equality of the ghost and 
quark field strength renormalisation factors cannot be expected. This relationship will be used as a measure of the 
fidelity of the model below. 

We remark that the form of Eqs. 



ff 



f 3 follows tradition in that vertex renormalisation factors are made explicit. 
However, this is in fact not necessary, or even desirable: once a vertex model has been adopted it can be regarded as 
renormalised (just like the quark dressing functions) and hence the renormalisation factors can simply be absorbed 
into the vertex factors V and W. Our prescription is equivalent to making this choice. 

Once Eqs. |f6| and[f7|are in place, renormalisation can be completed by performing a single subtraction at a scale 
fi to determine Z c (n, A) and Zf(/i,A). We remark that the subtracted gap equations are often renormalised at large 
Euclidean momentum where it is assumed that the perturbative form of the quark propagator is recovered. In general 
this is not true, and we choose to retain the constants -B(/x), A(p), and h(/i) as possible input parameters of the 
theory. 



C. Numerical Results and Comparison to the Lattice 

The primary ingredient in this investigation is the gluon dressing function, G. The unrenormalised propagator has 
been computed in Landau gauge in pure SU(3) gauge theory on large lattices [T4"j. with results shown in Fig. [3] 
A good fit to the lattice data is obtained with 



our) = 1^09 '-^-^f- 1 - 144 " GcV- 

p z — 1.144 



(18) 



The figure also displays an alternative fit of the form 
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FIG. 3: Lattice results fH] for the Landau gauge gluon propagator and model fits. Error bars are smaller than the points. 
Dashed line Eq. [T8j dotted line Eq. [19] 



which was specifically chosen to be similar to Eq. 18 in the 'Euclidean region' (p 2 E > 0) but to differ dramatically in 
the 'Minkowski region' (p 2 E < 0). Notice that the expected asymptotic behaviour G — > 1/p 2 + logs is not incorporated 
in the fits. This should not affect glueball properties or the long range structure of the propagators since these are 
dominated by the infrared regime. Although the lattice data carry units, the propagator has not been normalised 14 ; 
we therefore must allow for a multiplicative renormalisation factor in subsequent computations involving the gluon 
propagator. Finally, as indicated in the figure, these fits yield essentially identical results in the Eucli dean r egion, but 
differs crucially in the Minkowski region, p 2 < 0. The importance of this will be discussed in section MB 



The ghost propagator is computed next since it may be obtained self-consistently with an Ansatz gluon propagator. 
Even though the ghost lattice data are unrenormalised, we chose to renormalise with respect to them so that a 
reasonable com paris on can be made. In this case we set /j, = 3.57 GeV and h(p) — 1.27. Lattice data and numerical 
are displayed in Fig. |ij We chose to scale the gluon propagator so that G#(^ g ) = 1/Mg> thus 



solutions to Eq. 



13 



Gr(p; fig) 



gigt(g) 
Gi at (%)M§' 



(20) 



The coupling g 2 ([i g ) is then fit to the lattice result to obtain the final numerical prediction for the ghost propagator. 
The result for the case of the bare ghost-gluon vertex is shown as a dashed line in the figure. In this case the gluonic 
renormalisation scale was set to 1 GeV and the fit coupling was g(l) — 2.4. Notice that since the coupling and gluon 
propagator occur together in the gap equation, the scale dependence of the coupling is simply given by 



9 Og) = 9 (1) r , (21) 

where /j, g is measured in GeV. 

The quality of the fit is quite good over the entire range of the propagator, although it does worsen near the origin. 
An alternative fit using the vertex model of Eq. [6] (or Eq. 14 ) is shown as a dotted line in the figure. In this case 



the propagator saturates in the infrared and does not produce the sharp peak seen in the lattice data. We therefore 
prefer to use the bare ghost-gluon vertex in the following. 

Finally, we consider the quark dressing functions. We compare to the lattice computation of Ref. |15j where 
the mass function M = A/B and the wavefunction renormalisation function Z — 1/A were computed in quenched 
and unquenched Landau gauge QCD with 2+1 flavours. In this case the mass function does not depend on the 
renormalisation parameters and a direct comparison can be made. Alternatively, the authors of Ref. |15j chose to 
renormalise A by setting Z(p — 3 GeV) = 1, which we also employ. 
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FIG. 4: The renormalised ghost dressing function and lattice resultsJT?]. Error bars are smaller than the data points. Dashed 
line: bare vertex model; dotted line: full vertex model. 
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FIG. 5: The quark wavefunction renormalisation, Z, and quark running mass B, along with lattice results 15 . 

Results of the computation are shown in Fig. [5] The lattice data are unquenched results in the chiral limit. The fit 
coupling is g(l GeV) = 4.8 and the modified central Ball-Chiu vertex with n = 2 has been employed. We have found 
that fits with n = 1 give similar results but with larger couplings. Again, the agreement is satisfactory, although the 
wavefunction renormalisation deviates from the lattice data by 20% in the infrared region. 



III. GLUEBALL BETHE-SALPETER EQUATION 
A. Notation and General Expressions 



The Bethe-Salpeter equation is traditionally derived by considering four-point scattering near resonance locations 
via an equation that can be written schematically as 
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T 4 = V T R ^-T R + regular. (22) 

a s - m R 

The vertex factors are the Bethe-Salpeter amplitudes that describe the bound state resonance in question. In the case 
of gluon-gluon scattering the BS amplitude is given by 

T%(x,y,P) = (0\TA^x)A"(y)\R(P,X)) (23) 

where the ket labels the glueball in terms of its four-momentum and (possible) polarisation index. The Fourier 
transform of the BS amplitude is defined as 

r^( X ,y;P)=e-* p -^ J 0^ e ~«-^-v) r"" (k, P) (24) 

where k is the relative four-momentum describing the two-body state. It is convenient to introduce a modified BS 
amplitude, Xj via 

IV(fc, P) = D lta (k+)x a p(k+, k-)D pv (k-), (25) 

where 

k± = | ± k. (26) 

Bethe-Salpeter equations are typically simplified by projecting the BS amplitudes onto different Dirac or Lorentz 
structures. Thus, for example, the pseudoscalar quark amplitude can be expanded as[16 

X (k, P) = l5 {tE + jtF + PG + i\jt, P]H}. (27) 

We derive similar expansions by considering transversality and parity constraints on the gluonic Bethe-Salpeter 
amplitude. 

Transversality of the gauge field implies that fc+r^fc, P) = k^_T liV {k, P) = 0. If one considers the gauge £ = oo 
one obtains the transversality conditions 

k^u(k+,k.) = k v _ X »u{k+, = 0. (28) 

Application of the parity operator to the matrix element in Eq. [23] yields the following relationship 

y%(k+,kJ) = e(fx, v)-Jl R - xgfr+X) (29) 

where Tl R is the parity of the resonance R, k = (fco, —k), and e is a matrix of signs with values e(0, 0) = e(i,j) = 1 
and e(0,<) = e(i,0) = -1. 

The general Lorentz structure for a scalar glueball must be 

X""(fc+, fc_) - Aog^ + A x klk v + + A 2 k»k u _ + A 3 k»k v _ + A^_k\ + A 5 e" vaP k +a k^ . (30) 

Applying the transversality and parity constraints then yields the relationships 

X^ v (k + ,k.) = e^ a0 k a PpF(k,P) [0-+] (31) 

for the pseudoscalar glueball and 

X tlu {k+,k-)=A{k,P)A flu +B(k,P)B^ [0++] (32) 

for the scalar glueball. Here 



and 



(fc + ) 2 
k+ -k- 



fc+ • 



(35) 



(36) 



The Bethe-Salpeter equation considered here is displayed in Figs. |6j[7j and [8] Notice that quark and ghost Bethe- 
Salpeter amplitudes appear in the first equation (Fig. [6|. This may surprise the reader since no such amplitudes figure 
in the expansion of the scattering amplitude alluded to above (Eq. [22] ) . But the equations for the quark and ghost 
BS amplitudes can be substituted into the gluon equation to obtain the expression indicated diagrammatically in Fig. 
[9j This equation follows the usual BS conventions and indicates that the kernel is given by an infinite collection of 
diagrams that subsumes gluon exchange, the gluonic contact interaction, and dressed quark and ghost box diagrams. 

Although the resulting BS kernel is quite sophisticated it represents an extreme truncation of the diagrams that 
contribute to the gluonic bound state. For example, one can speculate that diagrams in which many gluons are 
emitted from the BS amplitude are important at long distances. Nevertheless, we shall press ahead to determine if 
this relatively simple model can yield a spectrum with any resemblance to lattice data. 




FIG. 6: Gluonic Bethe-Salpeter Equation. Dots represent full propagators and model vertices. Crossed diagrams are not shown. 
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FIG. 7: Ghost Bethe-Salpeter Equation. 



In the following the ghost BS amplitude will be denoted x, while the quark BS amplitude will be denoted 
the coupled Minkowski gluon-ghost-quark system in Landau gauge is given by the following equations. 



Then 



X fl u{k+,k-) = ig 2 N 
+ig 2 N 



d 4 q 
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FIG. 8: Quark Bethe-Salpeter Equation. 
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FIG. 9: Iterated Bethe-Salpeter Equation indicating the BS kernel in the pure gauge case. 
+i9 2 N J ^ iX {q+,q-)Q l ,Aq,k;P)H(q + )H(q_)H{Q) 

-y J -^ti[^S(q + )x(q+,q-)S(-q-) 7u S(Q)} h n (k + )h n (k_)A(Q, q+)A(Q, q_) 
+ crossed 



X (k+,k_) = ig 2 N I ( ' lq 



+ig 2 N 
= 9 2 Cf 



(2tt)* 

d 4 q 



+ig 2 C f 



(2-k) 
d A q 

d i q 



X (q+,q-)H(q, k; P) H(q+)H{q^)G{Q) 
4 X*p(q+, q-) B a ?(q, fc; P) G(q + )G(q_)H(Q) 



(37) 



(38) 



(27T) 



7a5'(-Q)7^(g + )G( g _)x^(g+,g-)/i n (g + )ft"(g-)A(Q,fc + )I(Q,fc_) 

7^(<?+M<Z+, «_)£(-«_)7„G(Q)iV(Q) /i 2 "(Q) A(k+, q+) A(k-,q-). (39) 



where Q = q - k and A(fc,p) = (^4(fc 2 ) + A(p 2 ))/2. 

The tensors in these expressions are as follows. The gluonic contact term is 



C a/? M » _ 2g^g af} - g av g^ - g^g 13 ". 
The three-gluon interaction with the vertex of Eq. [S] is given by 



T a p»v = P 11 (Q)V a/J ,y(q + ,k + )Vp v ^(q-,k-), 



(40) 



(41) 



where the projection tensor is 



Frr (Q) = g< 



Q~>Q~< 



(42) 
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and the three-gluon vertex tensor is 

V am {q +1 k + ) = A+(q + , k+, q+ - k+)g^{q + + k+y + A+(k + , q+ - k+,q + )g^(q+ - 2k+) a + 

+A + (q + - k + ,q + ,k+)g1 a (k + - 2q + r. (43) 

Here we reference the vertex factor, A+, of Eq. [8] or Eq. [9j 
The ghost-gluon coupling tensor is given by 

^-(«^-0)(^-'(£H)- 

Finally, the tensors appearing in the ghost BS amplitude equations are 



and 



B af} = ~k%ki 



( Hl+) 
\HQ) 



- 1 



(HQ) 
U(A-) 



HQ) 

Hq-) 



- 1 



, %+) ,\ (Hi-) , 
+ h(k+) ) ' ^ HQ) 



Hq-) 

h(k_) 



- 1 



(45) 



(46) 



B. Spin Zero Glueballs 



The spin zero (even charge conjugation) spectrum is obtained by substituting the pseudoscalar and scalar expressions 
for the gluonic Bethe-Salpeter amplitudes, Eqs. [31~]and[32] into Eqs. [37] -[39] Equations for the unknown amplitudes 
are obtained by projecting where appropriate and solving an eigenvalue problem of the type 

F(k;P) = X(P) J d\lC{k,q;P)F{q;P). (47) 

The eigenvectors and eigenvalues depend parametrically on the glueball momentum P| , which is varied to obtain the 
glueball mass at A(P| = — M 2 ) = 1. 

Recall that the BS formalism is developed with the Euclidean metric, which means that the mass condition just 
specified can only be obtained with imaginary Euclidean glueball momenta. This presents a potentially daunting 
obstacle since variables such as q± and k± become complex and hence dressing functions must be known over (a 
portion of) the complex plane. Note that this problem does not arise in the gap equations, where knowledge of the 
ghost and gluon Euclidean dressing functions is sufficient to obtain other two-point functions. Different approaches 
to this issue include diagonalising for Euclidean values of P 2 and extrapolating to the Minkowski region, or solving 
gap equations in the complex plane [IT]. As will be seen, extrapolation is not viable in this case and one must work 
with Minkowski values of the glueball momentum. Since we have chosen to use lattice data to model the gluon and 
ghost propagators, this implies that we must find a parameterisation of the data which reliably carries it into the 
region P| < 0. Of course this is impossible without additional information. We have thus made two models for the 
gluon propagator that agree in the Euclidean region but yield drastically different results in the Minkowski region 
(the dashed and dotted lines of Fig. [3]). 

We remark that the model of Eq. |18| postulates an entire gluon dressing function, which violates causality unless 
colour sources are strongly localised. But this is consistent with confinement, and hence entire propagators are 
sometimes used in hadronic model building|18j. 

The application of the formalism developed here to the pseudoscalar glueball reveals substantial simplification. 
Specifically, the antisymmetry of the indices in Eq. [31] guarantees that the contact contribution to the BS equation is 
zero. Similarly the coupling to the ghost BS amplitude is zero. This is anticipated because under parity and charge 
conjugation the ghost Bethe-Salpeter amplitude transforms as 



Xn{k+,k-) = Tl R XR(k+,k-) 



(48) 



and 



XR(k+,k^) =C RXR {k-,k + ), (49) 
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where Hr and Cr are the parity and charge of the state R. Thus spin zero states of negative parity cannot be made 
and even charge conjugation requires that the BS amplitude is even in k + ■ fc_ or k ■ P. 

Finally, although we consider the pure Yang-Mills theory in order to facilitate comparison to lattice computations of 
the glueball spectrum, the coupling to pseudoscalar quark-antiquark components is also suppressed. This is because 
the dominant term in the pion BS amplitude is known to be E (see Eq. [27|[5], thus ~%_ oc 75. Since four Dirac 
matrices are required to yield a nonzero value when traced against 75 , the coupling of the pseudoscalar glueball to the 
leading component of the pion is zero. Thus the structure of the pseudoscalar glueball is dominated by the three-gluon 
portion of the BS equations. The scalar glueball is more complicated: ghost and contact terms contribute and, as 
shown above, there are two gluonic components in the BS amplitude. 

The coupled BS equations have been solved using a number of strategies. The simplest approach assumes the 
'angle approximation' wherein all angular dependence of the kernel, )C, is neglected. Remarkable simplification ensues 
and the eigenvalue problem can be easily and accurately solved. The angle approximation then serves as a useful 
diagnostic for the full problem. The latter has been solved by making an expansion in Chebyschev polynomials of 
the second kind (described in the Appendix). As an alternative, we also discretised in the (y, k) plane and directly 
diagonalised. We found that the partial wave expansion converged quickly, requiring about the same number of waves 
as the order of the eigenvalue. Also, the angle approximation was very accurate for the ground state. 

The solution strategy followed that for the gap equations: the gluon and ghost lattice propagators were multiplica- 
tively renormalised at a scale /i g and the coupling g(fi g ) was determined by fitting to lattice gauge results. For the 
latter we relied on Ref. [T5] who determined 

M(0++) = 1730 ± 94 MeV M(0++') = 2670 ± 222 MeV (50) 

and 

M(0" + ) = 2590 ± 136 MeV M(0" + ') = 3640 ± 189 MeV. (51) 

Statistical and scale errors have been combined in quadrature. It is of note that the newer work of Ref. |19j does not 
quote excited state masses, so perhaps these should be taken with caution. 

Given that the lattice results are not renormalised, it is tempting to implement a Wilsonian renormalisation scheme 
by imposing a regularisation with scale A and determining g(A) (and m(A) if present) by fitting to observables such 
as glueball masses. Unfortunately we have found that the gap equations are numerically unstable when implementing 
this scheme, and it is necessary to render the equations finite by renormalising. 

A number of surprising features emerged from the solutions. First, the effect of the ghost channel is very small 
and it can safely be ignored. This is perhaps related to a second peculiarity we observed, namely the ground state 
mass depends very weakly on the strength of the contact term, even if that strength is varied over several orders of 
magnitude. This strange behaviour was confirmed analytically on very small discretisation grids. 




0.5 1 1.5 2 2.5 3 3.5 4 
IPI (GeV) 



FIG. 10: Glueball Bethe-Salpeter Eigenvalues. Lattice data [19] are represented by the horizontal bars. Solid lines: ground 
states; dotted lines: first excited states. (/1 = 1 GeV, g(fi) = 0.51, intermediate vertex model). 

The sensitivity of the gluon two-point function to assumptions motivates examining a variety of three-gluon vertex 
models. We have found that the ground state scalar glueball mass is largely independent of the assumed form of the 
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vertex, however, the pseudoscalar does depend on the vertex model. If one set the coupling with the scalar mass, 
then using the bare three-gluon vertex yields a pseudoscalar mass of 2270 MeV. Employing either the von Smekal et 
al. or Pennington- Wilson vertex model gives 2990 MeV. The best result is obtained if an intermediate three-gluon 
vertex is used where the ghost propagators are kept at their bare values in Eq. [8] Fitting the scalar glueball mass in 
this case gives g(l GeV) = 0.51, a pseudoscalar mass of 2750 MeV, and the results shown in Fig. 10 We find that 
the excited scalar glueball is approximately degenerate with the ground state pseudoscalar. This also held when the 
bare vertex model was used, and is consonant with lattice results. However, the excited pseudoscalar state is too light 
with respect to the lattice, indicating that the effective interaction in this channel is too soft. 



All of these results have been obtained with the entire gluon propagator model. The alternate fit (Eq. 19 1 does not 
yield physically acceptable results: the eigenvalues never cross unity for reasonable couplings and are not monotonic 
in the coupling strength - indicating that the solutions are anomalous. It is thus clear that the spectrum is sensitive 
to the form of the analytic continuation of the gluon propagator into the Minkowski region and to the vertex models 
employed in the kernel. 



IV. CONCLUSIONS 



We have attempted to find a consistent description of the ghost and quark propagators and the low lying scalar and 
pseudoscalar glueball spectrum with the Schwinger-Dyson/Bethe-Salpeter formalism and the lattice gluon propagator 
as input. While a reasonably good description of all of the quantities was obtained, the couplings required to obtain 
such descriptions were drastically different (<?(1) = 4.8 for the quark propagator, g(l) — 2.4 for the ghost propagator, 
and g(l) = 0.51 for the glueball spectrum). Furthermore, recall that gauge covariance and our renormalisation choices 
require that Z c = Zp. Direct computation of these renormalisation constants yields Z c — 1.37 while Zp = 0.89 at 
fi g = 1 GeV. This 35% discrepancy again points to quite large residual truncation-dependence. It thus appears that 
a consistent description of simple properties of QCD remains to be achieved. In particular, the effective strength of 
the gluonic interaction in the model glueball kernel is much too large. 

A variety of theoretical and practical issues bedevil this work. Amongst the former are the breaking of gauge 
invariance and even renormalisability that occur when truncations to Schwinger-Dyson equations are made. Finding 
truncation schemes that retain fundamental properties of the Lagrangian is clearly a desirable goal. On the practical 
side, even simple two-point functions are strongly dependent on model assumptions. As an illustration of the issues 
involved, the Kugo-Ojima confinement criterion specifies that the ghost dressing function h approaches infinity at the 
origin [15]. The lattice (and our computation) thus do not agree with this requirement, implying that something vital 
is missing in our understanding of the infrared properties of QCD. 

There are, of course, more direct sources of possible trouble. For example, it is clearly desirable to compute the 
gluon propagator so that the analytic continuation to the Minkowski region can be made with some reliability. Vertex 
models need to be vetted and improved as much as possible. Gauge invariance needs to be maintained to the level 
of accuracy desired in the numerical outcomes. And the effects of topological or many-gluon terms in the interaction 
kernel need to be investigated. It will be of interest to include quark effects in the computation. Lastly, an attempt 
at noncovariant gauges such as Coulomb or axial gauge should be illuminating [2D]. 
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Appendix: Partial Wave Decomposition of the Bethe Salpeter Equation 



The generic BS equation is 



F(k,k-P) = / d 4 qlC(q 2 ,k 2 ,P 2 ;x,y,w)F(q,q-P) 



(A.l) 



with 
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k = (cos 0\, sin 9\ cos 62, sin B\ sin O2 cos #3, sin B\ sin #2 sin #3 ) (A. 2) 

and 

<f = (cos 01, sin 0i cos 02, sin 0i sin 02 cos 03, sin 0i sin 02 sin ^3). (A. 3) 

We proceed by taking P in the glueball rest frame and denning y = cos#i, x — cos0i, and w — q ■ k. It is possible 
(and crucial) to simplify this by choosing k in the 1-2 plane. Then w simplifies to 



w = cos f/i cos 0i + sin 0i sin tlx cos 02 = 



+ \J\-x 2 ^/\-y 2 z (A.4) 



with z = cos 02- With this choice the 03 integral is trivial and the 82 and #3 dependence drops out. 

The partial wave decomposition is defined in terms of Chebyschev polynomials of the second type which are 
orthogonal over a measure appropriate for the four-dimensional integrals being considered here: 



1 / 7T 

dx\f\-x 2 U n (x)U m {x) = -5 nm . (A. 5) 
We define 

F(k,y) = ^2u n (y)F n (k). (A.6) 

n 

With this decomposition one obtains 

F n (k)=Aj2 I q 3 dqF m (q)IC nm (q 2 ,k 2 ,P 2 ) (A.7) 

with 

fC nm = / dxdydz\J\ - x 2 a/1 - y 2 U n {y)U m {x) K.(q 2 , fc 2 , P 2 ; x, y, w). (A.8) 
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